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Abstract — In this work, we derive a new upper bound on the 
termination time of the Hegselmann-Krause model for opinion 
dynamics. Using a novel method, we show that the termination 
rate of this dynamics happens no longer than 0(n 3 ) which 
improves the best known upper bound of 0(n 4 ) by a factor 
of n . 

Index Terms — Hegselmann Krause Dynamics, Opinion Dy- 
namics, Distributed Rendezvous, Consensus. 



I. Introduction 

One of the key challenges in modeling of social interac- 
tions in a complex multi-agent environment is the modeling 
of opinion dynamics and how agents influence each other's 
opinion and how new technology and ideas diffuse through 
such a network. 

One of the models that addresses such dynamics is 
Hegselmann-Krause dynamics which is introduced in [5]. 
Since, the introduction of the model many attempts have 
been performed to estimate the termination rate of such 
dynamics and still, up to now, the exact termination time 
in the case of scalar dynamics is not know. 

Because of the simple nature of the Hegselmann-Krause 
model, this model inspired other engineering applications, 
especially those from multi-agent systems to apply this 
dynamics to a given problem. As an example of such 
applications is distributed rendezvous problem in a robotic 
network such as a network of space shuttles. In this problem, 
one may want to gather a set of robots which lack a central 
coordination to a common place. One approach to handle 
this problem is to use the Hegselmann-Krause dynamics in 
such a network. An overview of this method can be found 
in [3]. 

A lower bound can be derived for this dynamics, that 
shows that this system needs at least 0(n) iterations before 
termination. The best known upper bound is developed in 
[9] where it is shown that the termination time happens in 
at most 0(n 4 ) iterations. This bounds relies on a quadratic 
time-varying Lyapunov function which is developed and 
studied in [11]. 

Conventionally, majority of the previous studies in the 
domain of distributed averaging, either rely on diameter- 
type Lyapunov function [4], [13], [14], [6], [8], or rely on a 
quadratic Lyapunov functions [2], [7], [10], [12], [9], [11]. 
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In this study, we chose a different path to derive the 
new bound. The path to derive our bound is as follows: 
we consider a linear function of opinions of a well-chosen 
subset of agents. Unlike the usual Lyapunov functions, this 
linear function is not necessarily decreasing, however, we 
can upper bound the number of times it can increase, as 
well as the amount by which it can increase each time. We 
show that except an 0(n) steps, this function decreases by 
at least O(-) in all the remaining times. The fact that the 
total increment of the function is bounded above by 0(n 2 ) 
will conclude our proof. 

The structure of this paper is as follows: in Section[II] we 
briefly discuss Hegselmann-Krause dynamics and formally 
state the problem, with the main result of this work. In 
Section [TTTJ we review some known results, and present and 
analyze the Lyapunov function, which is the core of this 
work. Finally in Section [IV] we conclude our discussion 
with some directions for future studies. 

II. Problem Statement 

Consider n agents on the real line. We denote the opinion 
of agent i at time t by xt(t) G M, and so x(0) = 
(xx{Q), £2(0), x n (0)) S M. n , represents the initial pro- 
file of the agents. The opinion profile of the agents evolves 
over time according to the following dynamics: 



Xi(t + 1) 



1 



for t > 1, 



where 



Afi(t) = {j e [n] : \ Xj (t) -Xi(t)\ < e} 



is the set of neighbors of agent i at time t, and we 
used the short-hand notation [n] = {l,2,...,n}. Note 
that the opinion dynamics is completely determined by the 
confidence values e, and the initial profile x(0). 

The interesting object in this work is the profile of the 
system at time t, i.e., x(t) = (x\(t), X2(t), . . . , x n (t)), and 
its dynamic over time. It is shown in [1] that this system 
converges to a steady-state x* in finite number of steps, i.e., 
there is a finite time instance T after which the profile does 
not evolve any more. The termination time of the dynamics 
is formally defined as 

T = min{f : x(t) = x(t'), W > t}. 



e 



e 




= {ari(t),...,Z| W(t) |(f)} x q {t) 



Fig. 1. Illustration of notation. 



Our goal in this work is to find an upper bound on 
T, To this end, we define a Lyapunov function for the 
dynamic, and study the variation of this function over time 
for the process. This allows us to bound the number of steps 
its takes for the dynamic to achieve its steady-state. The 
following theorem states the main result of this work. 

Theorem 1: For any scalar Hegselmann-Krause dynamics 
we have 

T <3n 3 + n. 

III. The Lyapunov Function 

We need to define a some notation and review some 
known results before presenting and analyzing the Lyapunov 
function. Fig.Q]is a pictorial representation of the following 
definitions. 

Without loss of generality we may assume that the agents 
are labelled so that their initial profile is in an increasing 
order, that is xi(0) < ara(O) < ■ • • < x n (0). It is easy to 
verify that this order is preserved over time, as precisely 
stated in the following lemma. 

Lemma 1: The order of the opinion of the agents is 
preserved over time, i.e., x\(t) < X2(t) < ■■■ < x n (t), 
for all time indices t. 

Definition 1: We denote by U(t) the set of agents with 
agreed opinion on the minimum value among all the agents 
at time t, i.e., 

U{t) = {i : Xi (t) < x 3 {t), Vj = l,...,n}, 
It is clear that once two agents share the same opinion at 
time t, they share the same neighborhood A/"(t) for the next 
time instance, and therefore they follow the same dynamics 
from time t onward. The following lemma presents this 
property of U{t). 

Lemma 2: The sequence of sets {U(t)} is increasing over 
time, i.e., 

W(0) C U{1) C U{2) C • • • C U{T). 

Let v(t) = \U(t) \ + 1, which implies x v u\ is agent with 
smallest index whose opinion is located at the second most 
left position on the real line at time t: 

v{t) = mm{j : Xj(t) < x^t), Vi G [n]\U(t)}. 



Next, we defineQ T(n) = U{n) U {x v{t ){t)}- This is the 
set of objects who define the Lyapunov function as follows. 

Definition 2: For the scalar Hegselmann-Krause dynam- 
ics, we define the Lyapunov function as 

£(*) = (x n {t) - Xi(t)) 

»eT(t) 

= |Z/(t)|(s n (t)-xi(t)) + (x n (t) -a;„ (t) (t)) (1) 
We can bound the value of C(t) at the initial and steady 
states. We also study the behavior of C(t) as dynamics 
evolves to find an upper bound for the termination time of 
the dynamics. 

In the following section, we first focus on dynamics with 
singular point steady-state, i.e., we assume that at the ter- 
mination time we have x\{T) — 22 (T) = • • • = x n {T). We 
show that the termination time for such dynamics satisfies 
the claim of Theorem 1 . Then we generalize this argument 
to cover systems arbitrary steady-state in Section IIII-BI 

A. Dynamics with Singular- Opinion Steady-State 

In the following we assume that the profile of the agents 
at the termination time is single point, and the system never 
splits into isolated sub-systems. First note that 

£(0) < i/(0)(x„(0) - n(0)) < v(0)ne. (2) 

We also point out that the value of the function at the steady- 
state is a finite and non-negative real number. If the system 
converges to a single opinion (as the focus of this section 
is), then £(T) = 0, but £(T) > in general. 

It is worth mentioning that C(t) inherently depends on 
longest distance between the agents in the system, as well as 
the number agents collapsed at the most left position. While 
the former decreases over time (Lemma[TJi, the latter (\U{t)\) 
is increasing. As a subsequence, the Lyapunov function 
is not a monotonic function over time. Fig. [2] shows an 
example of evolution of C(t). However, each increment in 
C(t) corresponds to merging of at least a agent to U[t), 
and therefore since this can only happen n times within the 
entire process, C(t) is decreasing unless for at most n time 

1 Note that in general there may be more than a single agent collapsed at 
the most second most left position. However, T(t) includes only the one 
with smallest index. 



instances. To make it more precise, 

1 = {0 < t < T : v(t + 1) > v(t)}, 
V = {0 < t < T : u(t + 1) = v(t)}. 

In the following we analyze the variation of C(t) for t £ I 
and t € T>, separately. For { £ I, although the Lyapunov 
function may increase, we show that the total amount of 
increment for all t € I cannot exceed a certain upper 
bound. On the other hand, for t £ T> the Lyapunov function 
decreases, and the amount of decrement is lower bounded. 
Hence, the number of time instances in T> is bounded from 
above. 

a) Variation of C(t) for t £ I: For t £ I, we have 
v(t + 1) > v{t) + 1. Let M(t) = Ni{t) \ T(t), m(t) = 
\M(t)\, and x(t) = ^£ teM(t) ar 4 (t).Itis clear that x t (t) > 
x v{t) if) for i G M(t), and hence x(t) > x v ^ (t). Therefore, 



x x {t + l) = 



E 



i£7Vi(*) • hi 



(*) 



> 



|M(t)| 

j^WM*) + ^( t )(t) + E tg M(t) ^0) 

|W(t)| + l + |M(t)| 
|W(t)|a:i(t)+a: K t)(t) 



(3) 



|W(*)I + 1 

which implies 

\U{t)\ Xl {t)+x v(t) {t)< V {t) Xl {t 

Therefore, the difference between C(t + 1) and C(t) for 
t &1 can be upper bounded by 

C(t + 1)-C(t) 

= v{t+l)x n [t+l)-\U{t + l)\x x {t+l)-x vm) {t+l) 
- v{t)x n (t) -\U{t)\ Xl {t) -x v(t) {t) 
v{t + \)x n {t+\)-v{t)x n {t) 

- v{t + l) - v(t) x^t + 1) 

- Xi(t + 1) - X v (t + t)(t + 1) 

u{t + l)x n (t + l)-u(t)x n (t + l) 

- u(t+l) - u(t) x t (t+ 1) 
u(t + 1) - v(t)\ \x„(t + 1) - xi(t+ 1) 
v(t+l)-v(t) 



< 



< 



(4) 



(5) 



< 



(6) 



where in (HJi we used the inequality in (fjjl, and (0 is due the 
facts that x n (t) > x n (t+ 1) and + 1) < x v u + u (i + 1). 
Hence, we have 



£[£(* + !)-£(*)]<£[»/(* + !)-«>(*) 



tei 



= ^(T) — i^(O) J ne = n-v(0) 

2 Recall that {u(t)} is in increasing sequence 



ne (7) 



\C{t) 



£(0) 





t 


*1 


«1 *3 < 


4 


*5 


*6 


T 



Fig. 2. Evolution of For any dynamic with singular-opinion steady- 
state, C(T) = 0. The function is always decreasing except at agents X = 
{t 1 ,t 2 ,...,t 6 }, with |Z| <n. 



Therefore, using (f2|i we get 

£(0)+5^[£(t+l) -£(*)] <n 2 e. 

£>) Variation of C(t) for t G P." In the following 
we study the evolution of C(t) for t G T>. First note the 
following lemma. 

Lemma 3: For any t £ T>, we have Af v u\(t) \Afi(t) = 0. 
Proo/: First, it is clear that M(f) C A/" v ( t )(t), other- 
wise the distance a;„( t )(i) — a;i(i) > e, and hence the entire 
system is split into two isolated systems U(t) and [n] \U(t), 
and cannot converge to a singular point. 

Now, assume the claim is not true, and Mi (t) — Af v (t) (t). 
In this case, we have 



+ 



ieAT„ (t) (t) 



|M(*)| 



(8) 



and therefore, G + which yields f(t+l) > 
Note that the latter is in contradiction with t £ T>. ■ 

Now, let q £ N v (t)(t)\Afi(t). Denote by d(t) = x v ^{t)- 
x\(t) < e the distance between the most left and second left 
set of collapsed agents. Let M(t) = Mi{t) \ T(t), m(t) = 
\M(t)\, and x{t) = ^ E ie M(t)^(*)- Note that X(t) > 
x v (t){t), because all the elements in M(t) have value which 
are not less than x v n\. Hence we have 



xi{t + l) 



EieM(i) 



Xi{t) 



> 



|M(t)| 

Eig^(t) x i(t) + x v{t)(t) + EigM(t) 

+ m(t) 

\U(t)\ Xl (t) + (xi(t) + d(t)) + m(t)X(t) 
v(t) + m(t) 

U{t)\ Xl (t) + (zip) + d(t)) + m{t){xi_{t) + d(t)) 

v(t) + m(t) 

, . m(t) + 1 ,,,, 
aci(t) + . 7 — ^—d(t). 



The same variation holds for all the agents in U{t), since 
they all share the same influencing neighbors, and hence, 



Hence, £(t) is decreased by at least 4- for each t £ T>. 
Therefore, 



x t {t + 1) = x 2 (t + 1) = • • • = x m) \(t + 1). 
On the other hand, 



T-l 



^(t)(i + 1) 



IMWI + l^ (t )W\M(i)| 



|Xi(t)| + l 

W |^(*)|a;i(t) + (m(t) + l)x u{t) {t) + (sift) + e) 

u(t)+m(t) + l 
= v{t){x v{t) (t) - d(t)) + (m(t) + l)x v(t) (t) + e 
+ m(i) + 1 
e - v(t)d(i) 



I " (i)(t) + Kt)+m(t) + l 



where (a) holds since the more agents exist in M v <t){t) \ 
Mi (t), the more moves :r„( t ) (i) towards the right. Moreover, 
in (6) we used the fact that x q (t) € M v (t){t) \ Mi(t), which 
implies x q (t) > xi(t) + e. 

Therefore, assuming v(t + 1) = v{t), we have 



£(t)-£(t + l) 

= [(Xn(t) - Xi{t)) - (x n (t + 1) - Xi (t + 1)) 



(<0 



2=1 
U{t) 



> ^2[Xi(t + l)-Xi(t)} 



, , , m(i) + 1 



v(t)d(t) 



u(t) + m(t) 



u(t) + m(t) + 1 



> 



v{t) + m(t) + 1 



(K*)-l)(m(t) + !)-!/(*) 
+ m(t) 



d(t) 



where (c) holds since x n (t) > x n (t + 1). Note that the 
second term above is strictly positive provided that m(t) > 
0. Therefore, for each t £ 2? with m(£) > 0, the Lyapunov 
function decreases by at least e/n. 

It remains to study the excluded case, where m(t) = 0. 
In this case the decrement in the Lyapunov function can be 
lower bounded by 



C{t) - £{t + l) > 
1 



v{t) + 1 



"(*) - 1 



it) + i 



v(t) v{t) + 1 

1 



d(t) 



(d) 
> 



/(t) + 1 



/(t) 



- 1 ( £ ) 1 

v{t)(v(t) + 1) £ ~ 3n £ 



where (<i) follows from the fact that d(t) < e, and the last 
inequality in (e) always holds for 3 < v(t) + 1 < n. 



0<£(T) = £(0) + £[£(t + l)-£(t)] 
t—o 

= £(0) + J] [£(i + 1) - £(i)] + £ [£(t + 1) - £(*)] 



< n 2 e 



which implies \V\ < 3n 3 . Hence, T = \T\+\V\ < n+3n 3 = 
0(n 3 ). Therefore, it takes at most 3n 3 + n steps to have 
£(T) > 0, which is the value of the function at the steady 
state for a singular-point steady-state system. 



B. Dynamics with Arbitrary Steady-State 

In this part we relax the assumption we made on the 
steady-state of the system in Section IIII-AI and show that 
the termination time of the system satisfies the claim on 
Theorem 1, regardless of the steady-state of the system. To 
this end, we split the entire process into several phases. The 
first phase includes the process from its beginning, until 
the most left agent achieves its steady-state. At this time, 
the system splits into two isolated sub-systems, namely the 
most left ones which remain constant for the rest of the 
process, and the remaining agents. The second phase starts 
at this time slot, and covers all the slots until the most left 
agent in the remaining set of agents gets to its steady-state. 
We can define the next phases similarly, and by adding up 
the duration of these phases we can bound the termination 
time of the entire process. 

Recall the argument in Section IIII-AI and recall that it 
remains unchanged for system with an arbitrary steady-state, 
except the claim on Lemma [3] in which we used the fact 
that Mi(t) C M u u\(t) holds over the entire process. This 
assumption violates if and only if x v ( t \{Ti) — xi(Ti) > e 
for some T\. In this case, we have Mi(Ti) = U(Ti), and the 
opinion of the agents in U(Ti) remains constant for the rest 
of the time, t > T\. In other words, the system splits into 
a time-invariant part U(T{) and the remaining [n] \U(Ti). 
Hence, we can proceed with the termination time of the 
remaining systems with n — \U(Ti)\ agents. 

As stated above, we denote by T\ the first time instance 
at which the agents in U(Ti) get split from the rest of the 
system. Now, consider the evolution of the profile from t = 
to t = Ti. During this interval, the entire argument in 
Section InLAl goes through. Define V x = X>n{0, 1, . . . , T x - 
1} and X x = in {0, 1, . . . ,T X - 1}. From © we can write 

J2 + 1) - £(*)] < H T i) - "(o)] ( 9 ) 

teii 

Also note that the Lyapunov function decreases by at least 



e/n for each t € X\. Therefore, 
< £(Ti) 

= £(0) + J2 W + 1) - + E ^ + !) - 
< v (Ti)ne-|X>!|^ (10) 
which implies £>i < 3^(Ti)n 2 , and hence 

Ti = |Ii | + |X>i| < i/(Ti) + 3^(Ti)n 2 . 

For time slots i > Ti, we redefine the Lyapanov function 
for the new agent set [n] \ U(Ti), with n - \U{T\)\ agents, 
and proceed until a new set of agents get isolated from the 
rest at some time slot T 2 . 

We can generalize this argument and define the following 
splitting-time set 

S = {Ti,T 2 , . . . ,T\ S \}, 

which is the set of boundary time indices of the phases 
defined above. We can further split the increasing and 
decreasing time sets as 

Z fe = Xu[T fc _i,Tfc), 
V k =VU[T k _ 1 ,T k ), 

where 

Pfc_i,Xfc) = {7fc_i,Tfc_i + 1, . . . ,7fc — 1}, 

and T = 0. Now, we repeat the above argument for phase 
the fc-th phase, and redefine the Lyapanov function for the 
system with remaining agents in this phase. In order to avoid 
confusion, we may use Uk{t) and C k (t) to denote the set 
of agents with opinion at the minimum value among the 
remaining agents in phase k, and the Lyapanov function in 
phase fc, respectively. We also define v k (t) — \Uk{t) \ + 1. 
Note that by definition, the system includes 

fc-i 

i=l 

agents in phase k, and this phase takes Tk—Tk-i time slots. 
Hence, we have 

< £k(T k ) = £k(T k -i) + + x ) - £fe W] 

tei,. 

+ + l)-£ fc (i)] 

tev k 

< v k (T k )n k e - \T> k \- — . 

3n k 

Therefore, 

\V k \ <2,v k {T k )n 2 k . 
This together with \I k \ < u k {T k ) implies 

T k - T fe _x < (3n 2 k + 1) v k {T k ) < (in 2 + l) v k {T k ). 



Finally, we can upper bound the termination time of the 
entire system by accumulating all such splitting times. 

\s\ \s\ 

fc=i fc=i 

where the last inequality follows from the fact that the 
total number of isolated points during the process does not 
exceed n. This concludes the desired result for dynamics 
with arbitrary steady-state. 

IV. Conclusion and Suggestions for Further 
Studies 

We studied the convergence rate of the Hegselmann- 
Krause model for opinion dynamics and using a novel 
Lyapunov-type function we proved a new upper bound of 
0(n 3 ) for its termination. From the practical perspective, 
this result is one step towards proving the scalablity of such 
a dynamics. 

An immediate extension to the current work is to use the 
technique introduced in this paper to prove a polynomial 
time upper bound for the termination time of the multi- 
dimensional Hegselmann-Krause dynamics which has been 
and is remained to be a very challenging problem. 
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